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TUNNELING FOR THE ROBIN LAPLACIAN 
IN SMOOTH PLANAR DOMAINS 

BERNARD HELFFER, AYMAN KACHMAR, AND NICOLAS RAYMOND 


Abstract. We study the low-lying eigenvalues of the semiclassical Robin Laplacian in a smooth 
planar domain with bounded boundary which is symmetric with respect to an axis. In the case 
when the curvature of the boundary of the domain attains its maximum at exactly two points 
away from the axis of symmetry, we establish an explicit asymptotic formula for the splitting 
of the first two eigenvalues. This is a rigorous derivation of the semiclassical tunneling effect 
induced by the domain’s geometry. Our approach is close to the Born-Oppenheimer one and 
yields, as a byproduct, a Weyl formula of independent interest. 

1. Introduction 

The spectral theory of the Robin Laplacian has attracted a lot of interest in the last years, 
especially in the strong coupling regime or, equivalently, in the semiclassical limit. Many authors 
have been interested in the asymptotic estimate of the bound states of this operator. The Robin 
Laplacian actually shares common features with the electro-magnetic Laplacian, the Dirichlet 
Laplacian on waveguides or 6 type perturbations of the Laplacian. These operators are often 
used to describe the physical properties of nanostructures (see for instance the review 0) 

In all these situations, numerous articles have revealed the role of the curvature in the cre¬ 
ation of eigenvalues or in the localization of the eigenfunctions. At some point, the case with 
Robin boundary conditions may also recall the boundary attraction that occurs for the magnetic 
Neumann Laplacian (and that is related to the surface superconductivity). At the scale of nanos¬ 
tructures the symmetries are known to induce a tunneling effect. This paper aims at quantifying 
this effect for bidimensional structures described by the Robin Laplacian on a smooth domain. 

1.1. Definition of the operator. Let 11 C M 2 be an open domain with boundary F = dQ. 
We will work under various assumptions on the domain FI. First, we consider the following two 
assumptions. 

Assumption 1.1. Fl is smooth with a boimded, regular boundary. 

As examples we can think of bounded domains (convex sets, annuli) or unbounded domains 
like the complementary of a bounded convex closed set. 

Assumption 1.2. The curvature k on the boundary T attains its maximum K max at a finite 
number N of points on T and these maxima are non degenerate. 

In the case when N = 2 in Assumption \1.2\ we will carry out a refined analysis valid under 
the following stronger (geometric) assumption: 

Assumption 1.3. 

i) Fl is symmetric with respect to the y-axis. 

ii) The curvature k on the boundary T attains its maximum at exactly two points ai and 02 
which are not on the symmetry axis and belong to the same connected component of the 
boundary. We write 

a 1 = (aqi, 01,2) £ T and 02 = (a.2,1,02,2) £ r , 
such that agi > 0 and 02,1 < 0. 

iii) The second derivative of the curvature (w.r.t. arc-length) at a\ and 02 is negative. 
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A simple example of a domain satisfying all the assumptions is the full ellipse 

f x 2 y 2 1 

< (x, y) : ^ ^ < 1 r ) witl1 0 <b < a. 

The two points in the boundary of maximal curvature are (±a, 0). The second example is the 
complement ary: 

f x 2 y 2 1 

<(x,y) : ^ + ^2 > 1 f > wit ^ 0 < a <b. 

The two points in the boundary of maximal curvature are (±a, 0). 

This paper is devoted to the semiclassical analysis of the operator 

U = -h 2 A, (1.1) 

with domain 

Dom(£/i) = {u, E H 2 (Q) : v ■ h%\7u — u = 0 on T} , (1.2) 

where v is the outward pointing normal and h > 0 is the semiclassical parameter. 

The associated quadratic form is given by 

VnGi^ 1 (fl), Qh(u) = f \KS7u\ 2 dx — hi f [u \ 2 ds(x), 

Jn Jr 

where ds is the standard surface measure on the boundary. 

Let (/i n (/i)) be the sequence of the Rayleigh quotients of the operator C\ l . It is known (see 
EH m M) that the bottom of the essential spectrum (if any) is non negative and this implies 
that, for all n E N, y n (h) belongs to the discrete spectrum as soon as h is small enough and that 
it is precisely the n-th eigenvalue of Ch counting multiplicities. 

The goal of this paper is to analyze the low-lying eigenvalues of the operator Ch in the semi¬ 
classical regime as h —> 0. The semiclassical analysis of the operator Ch naturally arises from 
the analysis of the Robin Laplacian with a large negative parameter a, 

^ — A Ra , Dom(—A r “)^ where Dom(—A^ Q ) = {u E H 2 (Q) : v ■ Vu + cm = 0 on T} , 

which has received a lot of attention (cf. [8, T7./U), jTJ) .TT|). The operator —A Ro! arises in several 
contexts, the long-time dynamics in a reaction-diffusion process [25] , and the critical temperature 
for enhanced surface superconductivity |12j . 

Putting a = — h~ 2 , we observe that a —> —00 as h —> 0 and the relation between the operators 
Ch and —A Ra is displayed as follows 

a(-A Ra ) = h~ 2 a(C h ). 

1.2. Known results. In this subsection, we recall the state of the art for this Robin problem, 
especially the spectral reduction of the operator Ch to an effective Hamiltonian on the boundary 
r. We will review an old result for the double well problem and apply it on the effective 
Hamiltonian. 


1.2.1. About the semiclassical Robin Laplacian. As a consequence of the results in [ I6l [T5l [32] 
we have the following theorem. 

Theorem 1.4. Under Assumptions QZ ^ and suppose that, among the maximal points of 
n, there are exactly M points ( a j)j£{i,...,M} where n" is maximal, then there exist a function 
h 1 —^ e(h) E (0, 00 ) such that 


lim e(h) = 0 , 
h->0+ 


and an interval 

I h = \-h- 


] -h-K max h 3/2 +'yh 7/4 -h 7/A e(h),-h-K miLX h 3/2 + ^h 7/4 + h 7/A e(h)[, 7 = \j R ^ , 

(1.3) 
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such that, for h small enough, 

a(C h ) n 4 = {hi( h),n 2 (h),. . .,HM(h)} , 

and 

m +1 (h) = ~h - K max /i 3/2 + 7 /i 7/4 + o(h 7/4 ), 7 = min f 37 , min ,/ ^ ( a j) 

Weaker versions of this result were obtained in [30] (and references therein, see also [27]). This 
result is related to B Theorem 1.1] in the magnetic case. 



1.2.2. About semiclassical tunneling on the circle. The aim of this article is to analyze the split¬ 
ting g .2 (h) — p>i(h) under the symmetry Assumption 1.3 (M = 2). We will see that the proof 
is easily reduced to the case when T has only one component (the one, by assumption unique, 
where n attains its maximum). As already observed in [15j the candidate for the splitting is 
obtained by considering the splitting for the operator 


= -h - Kma x h 2 + h 2 D 2 s + h -2 n(s), 
acting on the periodic functions in L 2 (M/(2L)Z), where 

l = H. 

2 ’ 


n = Krj 


— K . 


(1.4) 


and s the arc-length. Equivalently the operator can considered as the Schrodinger 

operator on the compact one dimensional manifold T. This is a double well problem which can 
be treated as a particular case of Helffer-Sjostrand [T8] with the effective semiclassical parameter 
being K := h*. 


Definition 1.5. We denote by /j® fF (/i) the j-th eigenvalue of Aif’ (counting multiplicities). 


Let us recall the splitting formula for the Schrodinger operator Aif rc ■= h 2 D 2 + t)(s) on 
the circle of length 2 L when V has two symmetric non degenerate wells at say s r and si with 
t)(s r ) = t)(s^) = 0 and D // (s r ) = V r, (se) > 0. We follow the exposition of |141 §4.5] (see also [35] 1 
but note that the formulas are established only for an example. In this paper, we will also use in 
many places the presentation of [2]. Because there are two geodesics between the two wells the 
discussion will depend on the comparison between the lengths of these two geodesics. For that 
purpose, let us introduce 


S = min (S 


u,S d ) , S u = f \/p(s)ds , S d = f \Jti(s)ds , (1.5) 

44,4 •/[«£,Sr] 

where \p, q] denotes the arc joining p and q in T counter-clockwise. 

The splitting formula for the operator A4f' rc is obtained by adding the “upper” and “lower” 
contributions and reads 

A2 (h) — Ai (h) = 4/Ptt- 2 7 2 ( AuV ^Me-ir 4. Ad \J t) {L)e~~h ^ + 0 {h 2e“"), ( 1 . 6 ) 

where 

(tZ) 7 (s) + 7 

[«r,0] \Msj 

(pl) / (s) -7 

hnM ZKs) 

7=(t» ,, (-r)/2)* = K(^)/2)* • 

Then, for the particular model A4^ ff , we easily notice that 

p.f(h)-pf(h) = h 3 H\ 2 (h)-\ 1 (h)), 




(1.7) 
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so that, under Assumption |1.3| we have 


t4 fF (h) - /if f (h) = 4/i 8 3 7r 2^2 f A u v /o(0)exp--^- + A d ^X>(L) exp- 

V h 4 /l4 / 


+ o 


U- + - 

ft* 4 exp - 


hi 


-7T ■ (1-8) 


Let us notice here that the complete proof of (1.6) provides a full asymptotic expansion and 
that the same holds for (1.8). Note that, if we assume that 0 is invariant under the symmetry 
exchanging the upper and lower parts, we have t)(0) = t>(L), S u = Sj and A u = A^. 


1.3. Statement of the main result. The main result of this paper is the following. 


Theorem 1.6. Under Assumptions]! . 1\ and\1.3\ we have 


t*2(h) - in(h) ~ u| {h)-nl (h) 

h-t 0 


(1.9) 


where ' [h) is defined in Definition 


formula in ( 1 . 8 ). 


1.5 


and where satisfies the asymptotic 


,eff 


The result in Theorem 1.6 shows a tunneling effect induced by the geometry of the domain 
(comparing with |19j . the boundary acts as the well and the points of maximal curvature as the 
mini-wells). This kind of reduction is also expected to be available for the magnetic Laplacian 
with a Neumann condition in smooth domains (see 13 Co]). However, magnetic fields induce a lot 
of additional difficulties especially in obtaining the optimal decay estimates of the eigenfunctions. 

Recently, magnetic WKB expansions are established in [Tj. Note that, in superconductivity, 
computing the splitting of the eigenvalues is useful to analyze the bifurcation from the normal 
state (cf. [TOj Lemma 13.5.4]). 

When the domain has corners and symmetries (e.g. the interior of an isosceles triangle), the 
tunneling effect is analyzed by Helffer-Pankrashkin in m • One difference between the setting of 
Theorem 1.6 and that in m appears in the spectral reduction to the reference problems. In na, 
the reference problem is a two-dimensional problem in an infinite sector which has an explicit 
groundstate. In this paper, the limiting reference problem is a direct sum of two one-dimensional 
operators. To prove Theorem |1.6[ we need to compare the eigenfunctions of the operator Ch 


with WKB approximate eigenfunctions (cf. Propositions |4.2| and 5.1). 

In higher dimensional domains, a spectral reduction, modulo 0(h 2 ), to an effective Hamil¬ 
tonian on the boundary is done in |32] (see also Section [7] where this reduction is explained). 
However, the analysis of the splitting as in Theorem |1.6| requires additional estimates since we 
want to control exponentially small error terms. 

In superconductivity, imposing a Robin condition, which is called in this context the de Gennes 
condition, models a superconductor surrounded by another normal/superconducting material (cf. 
|24l Theorem. 1.2]). In this context, we are naturally led to the analysis of the Robin Laplacian 
with a magnetic held where various regimes occur according to the comparison between the 
intensity of the magnetic held and the Robin parameter (cf. |21l 1221 20]). 


1.4. Organization of the paper and strategy of the proofs. Although it is easy to predict 
the statement in Theorem 1.6 once the effective Hamiltonian at the boundary is exhibited, the 
proof of the formula is much more technical. It will follow the steps outlined below: 


- In Section [2j we recall the known results related to the one dimensional situation. 

- In Section [3] we recall why the first eigenfunctions are localized, in the Agmon sense, near 
the boundary (the boundary is a well). As a consequence, we replace the initial problem by a 
problem in a thin tubular neighborhood of the boundary. Then the inhomogeneity of the new 
operator leads to a rescaling in the normal variable in the Born-Oppenheimer spirit and the 
introduction of the effective semiclassical parameter h = h 4 . 
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In Section [4j we analyze one mini-well problems (i.e. with one point of maximal curvature). 
Note that this terminology is the one of [19] where the problem was to analyze miniwells inside 
a degenerate well. We briefly recall the WKB constructions of |T6j in Subsection 4.2 Then, 


we establish optimal Agmon estimates in the tangential direction (see Subsection 4.3). 

In Section [5] we use the tangential estimates to prove that the first eigenfunctions are approx¬ 
imated (in the appropriate weighted space) by the WKB constructions. To do this, we are 
essentially led to use the same arguments as in dimension one (with respect to the tangential 
variable). Such estimates are closely related to the considerations of [28]| . 

In Section [6] we analyze the interaction between the mini-wells and establish Theorem 1.6 


Finally, in Section [7] independently of the tunneling problem, we derive a Weyl asymptotic 
formula for the counting function inspired by the considerations of [31] Chapter 13] and related 


to the effective Hamiltonian A4f". 


2. Robin Laplacians in one dimension 

Before starting the proof of our main results, it is convenient to introduce three reference 
lD-operators and to determine their spectra. These models naturally arise in our strategy of 
dimensional reduction and already appeared in [ 16] . 

2 . 1 . On a half line. As simplest model, we start with the operator, acting on T 2 (M_|_), defined 
by 

U 0 = -d 2 ( 2 . 1 ) 

with domain 

Dom(^o) = {u 6 H 2 (R + ) : u'( 0) = -u( 0)} . (2.2) 

Note that this operator is associated with the quadratic form 

r+oo 

Vo 3 u i—)• / \u' (t)\ 2 dr — \u(0)\ 2 , 

Jo 

with Vo = H 1 { 0, Too). 

The spectrum of this operator is {—1} U [0,oo). The eigenspace of the eigenvalue —1 is 
generated by the T 2 -normalized function 

uq(t) = a/ 2 exp(-r) . (2.3) 

We will also consider this operator in a bounded interval (0, T ) with T sufficiently large and 
Dirichlet condition at t = T. 

2.2. On an interval. Let us consider T > 1 and the self-adjoint operator acting on L 2 (0,T) 
and defined by 

nP = ~ d 2 , (2.4) 

with domain, 

Dom(?4 T} ) = {ii 6 H 2 (0,T) : u'(0) = -u(0) and u(T) = 0} . (2.5) 

The spectrum of the operator ?7 q ‘ is purely discrete and consists of a strictly increasing sequence 
of eigenvalues denoted by ^A n j . This operator is associated with the quadratic form 

3 u ha f \u'( t)\ 2 dr — \u(0)\ 2 , 

Jo 

with V ^ 1 ={«£ H 1 {0,T)\v(T) = 0}. 

The next lemma gives the localization of the two first eigenvalues Ai and A 2 (^ 0 ^) ^ or 

large values of T. 
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Lemma 2.1. As T -A +oo, there holds 

Ai(? 4 T} ) = -1 + 4(1 + o(l)) exp (-25) and X 2 ('hP)>0. (2.6) 

2.3. In a weighted space. Let 5 G M, T > 0 such that BT < Consider the self-adjoint 
operator, acting on T 2 ((0,T); (1 — 5r)dr) and defined by 

B { P = -(1 - Br)- 1 d T {l - Bt)8 t = -d 2 + 5(1 - BTpd T , (2.7) 

with domain 

Dom(? 4 T} ) = {«£ 5 2 (0,T) : u'(0) = -u(0) and u(T) = 0} . (2.8) 

The operator JJ-P is the Friedrichs extension in T 2 ((0,T);(1 — 5r)dr) associated with the 
quadratic form defined for u G vP = 5 1 ((0,T)) n {u(T) = 0}, by 

Qp( u ) = [ |w / ('r)| 2 (l — Bt) dr — |u(0)| 2 . 

Jo 

{T} 

The operator Ji B ‘ is with compact resolvent. The strictly increasing sequence of the eigenvalues 

r r n\ r rp'i f r j~'\ S 'TX 

of Wq 1 is denoted by (A n {n B J ) n eN*- It is easy to compare the spectra of J-l B 1 and Tif 1 as B 
goes to 0. 

Lemma 2.2. There exists To > 0 and C such that for all T > Tq, for all B G (—1/(3T), 1/(3T)) 
and n € N*, there holds, 

A n{U { p) - \ n (nP)\ < C|5|t( |A n (5f } )| + l) . 

Then we notice that, for all T > 0, the family (hlp^j is analytic for B small enough. In 

particular, its first eigenvalue Ai (j~Lp^ and the corresponding positive normalized eigenfunction 
{T} 

u B ‘ are analytic functions of B. 

Lemma 2.3. There exists Tq > 0 such that for all T > Tq, the functions (— 1/(3T), 1/(3T)) 9 
B h > Ai [hip^j and (—1/(3T), 1/(3T)) i-g uP are analytic. 

Proof. The family ('HP) does not fulfill the conditions for type ( B ) analytic op- 

V / B&(—Bo,Bo) 

erators in the sense of Kato since the parameter B appears in the definition of (the norm 
of) the ambient Hilbert space. Nevertheless, it becomes so after using the change of function 
u = (1 — Bt)~^u, since the new Hilbert space becomes L 2 ((0, T), dr), the form domain is still 
independent of the parameter and the expression of the operator depends on B analytically: 

UP = -(1 - BTpd T ( 1 - Br)d T ( 1 - Brp = -8 2 - ^ ^ )2 , (2.9) 

with the new Robin condition at 0 given by (0) = ( — 1 — y) u(0) and u(T) = 0. The price 
to pay is that the domain of the operator depends on B through the 5-dependent boundary 
condition. Note that the associated quadratic form is defined on 5 1 (0,T) by 

QpW = j Q IdrpdT- 4 (1 B Bt )2 dT ~ + f J 1^(°)1 2 • ( 2 - 10 ) 

□ 

r r p'\ 

The next proposition states a two-term asymptotic expansion of the eigenvalue Ai (fU B f ). 
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Proposition 2.4. There exists To > 0 and C > 0 such that for all T > To, for all B E 
(— 1/(3T), 1/(3T)) there holds, 




< CB 2 . 


{T} 

One will also need a decay estimate of u x B ' that is a classical consequence of Proposition 
of the fact that the Dirichlet problem on (0, T) is positive, and of an Agmon estimate. 


2.4 


Proposition 2.5. There exists To > 0, a > 0 and C > 0 such that for all T > To, for all 
B E (—1/(3T), 1/(3T)) there holds, 

II „ctr {7 1 } || . < O 

B L 2 ((0,T);(1—Br)dr) “ 

Remark 2.6. We will apply the results of this section with T = Dh~ r , r E (0, ^), B = and 
h E (0, ho) for ho small enough. 


3. Reduction to a tubular neighborhood of the boundary 

3.1. Agmon estimates. As proved in [IB], the eigenfunctions of the initial operator Ch are 
localized near the boundary and this localization is quantified by the following theorem: 

Theorem 3.1. Let eo E (0,1) and a E (CfyTo). There exist constants C > 0 and ho E (0,1) 
such that, for h E (0,ho), if Uh is a normalized eigenfunction of Ch with eigenvalue p < — eoh, 
then, 

[ (\u h (x)\ 2 + h\Vu h (x)\ 2 ) exp f - n dlst ( J 1 r ) \ dx < C. 

Jn \ h 2 J 

Hence, this theorem is a quantitative version of the statement that the boundary is a well (in 
analogy with the Schrodinger model in [IS]) as h —> 0. 


3.2. Spectral reduction. We can explicitly derive a reduction near each component of the 
boundary. From now on we assume for simplification that the boundary is connected. Given 
5 E (0, do) (with 5q > 0 small enough), we introduce the h-neighborhood of the boundary 

Vj = (x E h : dist(x', T) < h} , (3-1) 


and the quadratic form, defined on the variational space 

V$ = {u E H l (Vs) : u(x) = 0 , for all xEH such that dist(x, T) = 5} , 
by the formula 

Vu E Vs , Q { h 5} (u) = f |hV u\ 2 dx — h 2 f \u\ 2 ds(x). 

Jv s Jr 

Remark 3.2. In the following we will be led to take 6 = DhP with p E (0, . We will choose 

- either p < \ with D = 1, 

- or p = -^ and D > S where S is defined in (1.5) in order that the error term in (3.2) is smaller 
than the tunneling effect, we want to measure. 

Let us denote by p,j?\h) the n-th eigenvalue of the corresponding operator ■ It is then 
standard (cf. [18] ) to deduce from the Agmon estimates in Theorem 3.1 the following proposition. 

Proposition 3.3. Let eo E (0,1) and a E (0, yTo ). There exist constants C > 0, ho E (0,1) such 
that, for all h E (0, ho), 5 E (0, <5q), n > 1 such that p, n (h) < —eoh, 


Pn{h) < pffl{h) < p n {h) + Cexp aSh 2 ) . 


(3.2) 
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FIGURE 1. Illustration of the boundary coordinates (s,t) for the point x. Note 
that for t = 0 and s —> —L, the point x approaches the boundary point s = L. 


3.2 


3.3. Boundary coordinates. Thanks to 
. with the choice of 5 made in Remark 
near T we will use the canonical tubular coordinates ( s , t) where s is the arc-length and t the 
distance to the boundary. We recall some elementary properties of these coordinates. Let 


reposition 3.3, we can now work with the operator 
Since the functions of its domain are supported 


(M/2LZ) 3 s A M(s) £ T (3.3) 

be a parametrization of T (and thus we will always work with 2L-periodic functions sometimes 
restricted to the interval (—L,L]). The unit tangent vector of T at the point M(s) of the 
boundary is given by 

T(a) := M\s ). 

We define the curvature k(s) by the following identity 

T’(s) = n(s)v(s), 

where z/(s) is the unit vector, normal to the boundary, pointing outward at the point M(s). We 
choose the orientation of the parametrization M to be counter-clockwise, so 


det(T(s), i'(s)) = 1, Vs E (M/2LZ). 


We introduce the change of coordinates 

<I> : (M/2LZ) x (0,to) 3 (s,t) *->■ x = M(s) — tv(s) E V, 5 0 . (3-4) 

The determinant of the Jacobian of is given by 


a(s, t) = 1 — tn(s). 


(3.5) 


In the case of symmetry, we choose as origin of the parametrization the point po defined as 
follows 

{Po = (xo,yo),Po = (x 0 ,y 0 )} = T n {x = 0} and y 0 > y 0 , 
i.e. we suppose that s(po) = 0 and s(p o) = L. This is illustrated in Figure 1. 

3.4. The operator in a tubular neighborhood. We can now express the operator in these 
new coordinates. To indicate that we work in the coordinates (s, t), we put tildes on the functions. 
For all u E L 2 (Vs Q ), we define the pull-back function 

u(s, t) := u(<I , (s, t)). (3-6) 

For all u E H 1 (Vs 0 ), we have 



\u\ 2 dx = 


u(-s, t)| 2 (1 — tn(s)) dsdt , 


(3.7) 
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f \Vu\ 2 dx = f (1 — tn(s)) 2 \d s u\ 2 + \dtu\ 2 (1 — tn(s)) dsdt. 

Jv So J 1 

The operator c\^' s is expressed in (s, t) coordinates as 

c\^ = — h 2 a~ l d s (a~ l d s ) — h 2 a~ l dt{adt ), 
acting on L 2 {adsdt). In these coordinates, the Robin condition becomes 

2 3 

h dtu = —h^u on t = 0 . 


(3.8) 


We introduce, for 6 G (0, <5o), 

V <5 = {(s,t) : s G] — L, L\ and 0 < t < 5} , 

V s = {ueH 1 (V s ) : u(s,S) = 0 }, 

Vs = {u € H 2 (Vs) n Vs : d t u(s,0) = -h~^u(s,0)} , 

Q}P(u) = [a~ 2 \hd s u\ 2 + \hdtu\ 2 ^j a dsdt — hT J |rt(s, 0)| 2 ds , 

c\P = — h 2 a~ 1 d s (a~ 1 d s ) — h 2 a~ l dt(adt) ■ 


We now take 

5 = Dh p , (3.10) 

and write simply Ch for c\^ ' r . The operator Ch with domain V is the self-adjoint operator defined 
via the closed quadratic form V p 3 u Qh(u ) by Friedrich’s theorem. 


3.5. The rescaled operator. In order to perform the analysis and to compare with existing 
strategies, it will be convenient to work with a rescaled version of Ch- We introduce the rescaling 

(cr,r) = ( s,h~h ), 

the new semiclassical parameter h = h* and the new weight 

a(<r, r) = 1 — h^rn(a). (3-11) 


We consider rather the operator 


C h = h 1 C h 


acting on L 2 (ddadT) and expressed in the coordinates (<r, r). As in (3.9), we let 


Vt = {(<j, r) : a e] — L, L] and 0 < r < T} , 

V T = {m£ H l (yr) : u{a,T) = 0} , 

V T = {u € H 2 {Vr) n V T : d T u{a, 0) = -u{a, 0)} , 

Qh(u) = (d~ 2 h A \d a u\ 2 + \d T u\ 2 ^JddadT — J \u{cr, 0)| 2 da , 

CJ = —h 4 a~ 1 d C7 a~ 1 d CT —Ti~ 1 d T ad T . 

Remark 3.4. We then specify the analysis for 

T = h~U = Dh p "5 


(3.12) 


(3.13) 


and omit the reference to T. 
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4. Simple mini-well 


This section is devoted to the analysis of the eigenfunctions when the curvature has a unique 
non degenerate maximum (i.e. Assumptions 0 and |1.2| with M = 1). We will investigate both 


the WKB constructions and the accurate approximation of the eigenfunctions in such a situation. 
For that purpose, we will constantly work with the operator defined in the sequel. 

4.1. Definition of the simple mini-well operator. Let ui be an (open) interval in the circle 
of length 2 L identified with the interval (— L,L\. We can view w as a (curved) segment in the 


boundary of Q by means of the parametrization in (3.3). The operator h = C u is defined as 
follows. We assume that u contains a unique point s u of maximum curvature (i.e. n(s u ) = K max ) 
that is non degenerate. The form domain V u and the domain of this operator are defined as 
follows, 

i=wx (0, T), 

Vcj = {u € H l (Vj) : u = 0 on r = T and du x (0,T)} , (4-1) 

= {u £ H 2 (Vuj) fl Vu : d T u = — u on r = 0} . 

The operator is the self-adjoint operator on L 2 ( V u ; a da dr) with domain V w and 

= — h A a~ 1 (d rT a~ 1 )d C7 — a~ 1 (d r a)d T . (4.2) 

We denote by its lowest eigenvalue. 

Definition 4.1. The corresponding positive and L 2 -normalized eigenfunction is denoted by 4>h,uj- 

Let be the second eigenvalue of the operator C u . The analysis in [46] yields that, for 

h small, /-iuj(h) is a simple eigenvalue and 

p. 2 ,u(h) ~ h uj[K) = 3y hJ^ + ^/ 4 o(l) as K -> 0+ , (4.3) 


where 7 = \J ~ K ' . 

4.2. Reminder of the WKB constructions. 

4.2.1. Statements. In this section, we recall the WKB construction of |16| in the spirit of the 
paper by Bonnaillie-NoeKHerau-Raymond [I] (see also the classical references about the Born- 
Oppenheimer approximation [3] f26. 28]). 

Proposition 4.2. There exists a sequence of smooth functions (Oj) such that the following holds. 
We consider the formal series (or a smooth realization constructed by a Borel procedure) 

Vh,u( a ’ T ) ~ aj(cr,T) , (4.4) 

i>o 

where 

i) is the Agmon distance to the well at a = s u of the effective potential 


v(a) = K max - k ( ct ) 


and defined by the formula 


•Mcx) = / y/ti(a)da, 

J[su,c) 

ii) ao is in the form ao(cr, r) = £o,uj(c)uo(t) where 

uq(t ) = V2e~ T , 

and 

" a - 7 


£o,uj( cr ) = Co(o-) = (-) 4 exp 




da 
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is the solution of the transport equation of the effective Hamiltonian 

+ <9a($Uo) = , with 7 = \J K ^ . 

iii) For j > l, aj(a,r) is a linear combination of functions of the form 

fj,k{v)gj,k{T ), fj,k G C°° (ca) and g jjk G <S(M +). 

iv) The formal series 'I ' k w satisfies 

e*“ /H (X - /i) , = 0(h°°), 

where f is an asymptotic series in the form 

F - -1 - Kma x h 2 + 7 h 3 + ^2 Fjh j . (4.5) 

3> 4 

This series is the Taylor series of the first eigenvalue g u (h). 


In the previous proposition, we have used the following notation. 

Notation 4.3. We write a(a,T;h) ~ ^2j>o a j( a ' T )^ when for all J > 0, a 6 N 2 and all 

compact K C oo X M+, there exist hj^ a ,K > 0 and Cj, a j{ > 0 such that for all h G (0,hj, a ,K), we 
have, on I\, 


d a (a(cr,T;h) - ^ 


a,-(a, r 


i=o 


< Cj,a,i ^ J+1 


We also write a = 0(h°°) when all the coefficients in the series are zero. 


Remark 4.4. In the sequel, it will be convenient to work with a truncated version of''Let 
6 j be an open interval such that s u 6 w C ii C w and 

T ) = Xu>(<T)x(T- 1 T)Vh,u(<r , T ) j ( 4 - 6 ) 

where 

i) x is a smooth function cut-off function with compact support being 1 near 0; 

ii) Xu> G Cf°(oo) is a smooth cut of function satisfying 0 < < 1 and Xu = 1 on ui. 

The truncated function ifn,uj satisfies 

e*“ /h (Z u - f) ^h,u,u = 0(h°°) , in L \%), 


V j G {1, 2} , e*“/ h df (£„ - f) ^n,u,u = 0(h°°) , in L\%). 

In the sequel we will use that oj and ui can be chosen as large as we want, as soon as ui only 
contains one mini-well and ui satisfies the above condition. 


4.2.2. Proof. Let us just explain the main steps in the proof of Proposition 4.2 
formal Taylor expansion, we find the following expansion of the operator 


Thanks to a 


~ — d 2 — h^d 2 + 2h 6 Tn(a)dl + h 2 n(a)d T + h e Tn'(a)d a 


h 2j+2 T j (n{a)) j+1 d T - k'(<t) ^^ +6 d i r J (K((j))^ (T . 
j= 1 3=1 3=1 

We introduce the (formal) conjugate operator 


C h := exp 


d(o) 

h 


L h exp - 


d(a) 


h 
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and write 

- m) ~ 0 . 

We (formally) expand the operator £][ as follows 

OO 

1=0 

with in particular 

Qo = ~d 2 r , 

Qt = o, 

Qi = K(a)d T - d'(a] 2 , 

Qi = 2i9'(a)d a + r(a), 

Qi = -d 2 + C 3 T 3 k((T) 3 + T 3 (K(fj)) 4 9 r . 


(4.7) 


(4.8) 


We then rearrange all the terms in (4.7) in the form of power series in h 
and fi£ by expressing the cancellation of each term of the formal series, 
coefficient of H° yields the equation, 


and select $, a^(a,r) 
The vanishing of the 


(Qo ~ Mo)ao(o-,r) = 0. 


We have Q$ = Id <S> Ho on L 2 (M ct x R + , r ). This leads us naturally (considering the operator Ho 
introduced in ( |2.1| )) to the choice 

Mo =-1 and oq(<t,t) = £q(c)uq{t) . 


Since Qf = 0, the vanishing of the coefficient of h 1 in (4.8) yields 

(Qo ~ Mo)ai((r, r) - ma\(o, r) = 0 . 


This leads us to the natural choice fj i = 0 and 


ai(cr,r) = £i(cr)« 0 (r). 
We look at the coefficient of h 2 and obtain 


(Qo ~ Mo)fl2 + (Q2 ~ M2)ao — 0 . 

Remembering that /xq = —1 and ao(a,r ) = €o((?)uo(t), we get 

(Qo + l ) a 2 = -uo(t)(k((j) - d’(a ) 2 - n 2)Co(o") • 

By using the Fredholm condition with respect to r, we get the eikonal equation 

— n(a) — $ (a ) 2 — ^2 = 0 . (4-9) 

Consequently, we take /i 2 = —k(soj) and get d'(s w ) = 0 and we consider the solution such that 
^"(soj) > 0. This gives 


and 


*■<*.)= \J- K y . 

(4.10) 

/ V K max - n(a) da = $ u (a ), 

(4.11) 


where [s^, a\ is the segment joining s u and a counter-clockwise (the integral may also be under¬ 
stood as the Lebesgue integral on a measurable set, independently from the representation of 
the set). 
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We deduce that 02 is in the form 

02(cr, t) = &((t)uo(t) ■ 

Now we look at the coefficient of h 3 in (4.8). This yields 

(Qo + 1)03 + (Q2 ~ ^2)01 + (Q3 — 173)0.0 = 0 ■ 


Using (4.9), we see that the term (Q'] — ^ 2)01 vanishes and thus we get 

(Qo + 1)03 = ~(Qi ~ ^3)00 • 

For each fixed a, the Fredholm condition implies that 

((Q 3 ~ M3)ao, uo)l 2 (r+, t ) = 0 > 

that is 

W(v)€o(v) + (^"( 0 ) - M3)Co(o) = 0. (4.12) 

Since we look for smooth solutions at s u and for the smallest possible g. 3 , the linearization at 


a = s u , leads to 


h3 = 


~h"(Suj) 


= 7 . 


We can determine fo by solving ( |4. 12 ) in a neighborhood of er = and hnd 

1?" -7 


1 

7\4 


£ 0 ( 0 ) = - exp - 


7 r 


W 


-da 


where the constant is chosen to get a L 2 -normalized quasimode (modulo K). Then, we are led to 
choose 

03(0", t ) = 6(0')« o ( t ) ■ 

This construction may be continued at any order. 

4.3. Tangential Agmon’s estimates. In this subsection, we derive Agmon’s estimates for the 
eigenfunctions of the operator C u with domain and form domain introduced in (4.1). 
Let us start with the following elementary lemma that is related to the Born-Oppenheimer 
approximation. 

Lemma 4.5. There exist constants C > 0 and ho G (0,1) such that, for all h G (0, ho) and 
u G V u , 

Qu)(v) > / aT 2 h 4 |<9 t7 ’u| 2 adadr + / (—1 — n max h 2 + h 2 o(cr) — Ch 4 ) |u| 2 adadr . 

JVu, Jv„ 

Proof. Using ( 3.13[ >, we have 

r L r T 


Qw(u) = 


l-L JO 
,{T} u 


a z h A \d a u\ z + \d T u\ z )a drda — j \u(a, 0)| 2 da 


(4.13) 


Recall the operator in TF B s in (2.7). By a simple scaling argument and the min-max principle, 
we have 

f T f T 

/ \d T u\ 2 adr — \u(a,Q)\ 2 > \\{TL B ) / |u| 2 adr, T = h p ~ 2 , B = h* n(a). (4.14) 

Jo Jo 

Thanks to Proposition |2.4| we deduce the lower bound since D = K max — k. □ 


From Lemma 4.5 we may deduce some accurate tangential Agmon estimates satisfied by 
We will often use the following notation. 

Notation 4.6. For g G (0, L), we let 

B u (g) = (-Q+ Scj,q + s u ) and B w (g) = B u (g) X (0, T). 

Let us first state a proposition that will be convenient in the sequel. 
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Proposition 4.7. Suppose that T = Dh 4 and D > S. Let <t> be a Lipschitzian function that is 
a subsolution of the eikonal equation: 


D(cr) — |4 >/ ((t)| 2 > 0, VerGu;, 

and let us assume that there exist a non decreasing function M + 3 R 1-3 M ( R ) G R + tending to 
+00 as R — > + 00 , a positive constant ho such that, for all h £ (0, ho), and R > 0, 


O(o-) - |4 ,/ (cr )| 2 > M(R)h , 
|$(a)| < AL(R)h , 


V<7 £wfl CB u (Rh2 ), 
V<r G B u (Rti*). 


Then, there exist Ro,C > 0 and ho G (0, ho) such that the following holds. For all R > Rq, 
Co £ (0, 9 ) ), h G (0, ho)> -2 £ [~1 — ^max^j —1 — K max h 2 + Coh 3 ], U G T^uj, 


fi 3 ||e*Nl L 2(v„) < C'||e*/ ft (4 - z) u|| La( p u) + Ch 3 


\'n\\ \ 

1 "L*(y„pB w (im)) ’ 


and 


Z t 2 ^, , + Ch 3 \\u \\ 2 ^ . 

i2 (v u ) N L 2 (v u nB w (m2)) 


h 4 ||d CT (e #/ V)|| 2 2( ^ ) < Ch 3 \\e® /h (Cu - z^u 
Proof. By the usual Agmon formula, we get 

(Aj u,e 2 ®/ h u) = Qu,(e^^ h u) — h 2 f a _ 2 | < h / | 2 e 2 $ / r '|u | 2 ddadr . 

Jv„, 


(4.15) 

(4.16) 


By Lemma 4.5, we deduce 


(£u)U, e 2 ®/ h u) > f a 2 h 4 \d a (e®C l u )\ 2 ddadr 

JVu, 

+ f (—1 — K max h 2 + h 2 (p — a -2 ^'! 2 ) — Ch 4 ) \e^C u \ 2 ddadr . 

JVu, 

Note that, for all ( a,r) G V u , \h 2 K(a)r\ < Dh 4 p \K\oo. Thus there exists D, ho > 0 (depending 
only on p, D and |k|oo) such that, for all h G (0, ho) and (a, r) G V u , 


i ~ 2 - l| < Dh Ap . 


This leads to choose 


1 


P = 


and to the lower bound 

((Cu — z)u, e 2 ®l h u) > f d~ 2 h 4 \d a (e‘ !> Cu)\ 2 ddadr 


+ J 1 — n ma , x h 2 + h 2 (p — I4*'! 2 ) — Ch 3 — z^j |e $ ^ ri tt | 2 ddadr . 


for some given constant C > 0 independent of R. 
Using the assumption on z, we deduce that 


i- 2 h A \d a {e*/ n u )\ 2 adadr + (h 2 ( o - |T'| 2 ) - Ch 3 - ^p-h 3 ^J \e® /h u \ 2 adadr 


2 
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Now we use the assumption on the function $ and obtain 

( 


a 2 ft i \d rT (e®/ h 'u)\ 2 adcrdT+( —f-—^ — cA f h 3 \e^^ h u\ 2 adadT 

2 J JV W \1B U {RK L / 2 ) 


< ||e*/ fi (4, - z)u\\\\e* /n u\\ + Ch 3 \e^/ h u\ 2 adadr . 

Jv^rB^iRh 1 / 2 ) 

We choose Rq > 0 such that —C> 0. For all R > Rq, we have — C> — (7 > 0. 

Thus, by the Cauchy-Schwarz inequality and the assumption on the function 4>, 

J ^ d~ 2 h 4 \da(e^^ h u)\ 2 adadr + ^ ^ ^ ft 3 ||e $ /\t || 2 


< C/T 3 ||e $/s (iL--M | 2 + C/i 3 |M | 2 ^ . i 

i 2 (v u nfi w (Rfi2)) 


(4.17) 


From (4.17), we get 


/ MQRq) 


\ 3 

-C J /i§||e*/ ,i u|| < (7/i _ 5 || e $ / s (£ CJ — z )u\\ + Chi \ 


u\\ l 

N L 2 (v w nS w (.Rft 2 )) 


and we deduce (4.15). The estimate (4.16) directly comes from (4.17). 


□ 


Remark 4.8. If we apply Proposition f.l to the eigenpair (<fu},h, Pu{h)), wo get 

< C ||0£j,n|| L 2(y w ), 

(as soon as M is large enough, to insure that p u (h) belongs to the energy window). Note also 
that these estimates are weighted estimates in H)(uo, L 2 (0,T)). 

Let us gather some possible choices for <F in the following proposition (see |4, Chapter 6 ] or 
12 Proposition 2.4 and Lemma 2.5] for a detailed proof). 


Proposition 4.9. Let Co > 0 such that 

*>{&) > c 0 (a - s^) 2 and ^(cr) > c 0 (u - s u ) 2 , Vcr £ uo . 
Possible choices of 4? satisfying the assumptions in Proposition \f. 7| are: 

(a) for a £ ( 0 , 1 ), the rough weight 

<F = Vl - a 

with R > 0 and M = c^aR 2 ; 

(b) for N £ N* and h £ (0,1), the accurate weight 


(4.18) 


®u,N,n = ®u> - Nh In ( max ( —, N 




(4.19) 


with R = \ and M = N inf w 1 

(c) for a G (0,1), uo as above, TV G N* and h G (0, 1), the intermediate weight 

®u>,ui,N,h(p) = min \ ^uj,N,h(^), Vl~a inf 4> w (f) + / \/o(d) da > , 

[ suppXi, y J[cr,t] J J 

with R = and M = N min ^a,inf w where we recall that supported near du ; 

if Co is large enough. 

Moreover, the weight $ u ,u,N,h satisfies the following. Let K be a compact with I\ C [xdi = !}• 
For all N £ N*, there exists ao such that for all 0 < a < ao, there exist ho > 0 and R > 0 such 
that, for all h £ (0,ho), we have 

(i) ®u,&,N,h < on oo, 

(h) &u,ca,N,h = $cj,N, h on K, 
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(iii) $ 0 ) ,u,N,h = \/l - onsuppyL (c/. Remark 4 . 4 ). 

Remark 4.10. We now assume p = \ and D > S. 

5. WKB APPROXIMATION IN THE SIMPLE MINI-WELL CASE 

The aim of this section is to approximate the true eigenfunction q by the WKB function 
T/j u defined in (4.4), or equivalently by defined in (4.6). 

5.1. Main result. Let us introduce the orthogonal projection on the space spanned by 

— {i/j, (frhjUj ■ 

Recall that x £ C(?°([0,oo)) denotes a cutt-off function which is equal to 1 near 0. 
Proposition 5.1. Let K be a compact set in co. There holds 

(5.1) 

^ /h dA^ - nJ = 0(h °°), (5.2) 

in C(K; L 2 (0, T)) and where we have let = x(T 1 t )^h,oj- 

We split the proof into four steps inspired by the presentation of [2], We choose 6j so that 
K C {xuj = 1}. Then, it is sufficient (see Remark |4.4[ ) to prove that 

e*" A (^-n w VvO =o{h°°), 

e^ h dMh,u, -B^m) =0(h°°), 

in C(K; L 2 (0, T)). 

5.2. Estimating the L 2 -norm. Recall the definition of the domain in (4.1). Let 

u = il>n,u ~ n ■ 


(5.3) 

(5.4) 


(5.5) 


Since u is orthogonal to the eigenfunction 4>n,oj, then by the min-max principle, 


- Pu,(h))\\u\ 


L 2 (Vu) 


< 


2 

T 2 (W) 


( R- jj /I jj ( h ) ) Lh.Cj 


L 2 (Vui) 


Using the estimate of the P 2 ,u(R) ~ hi ( |4.3| ), the expansion of p w {h) in ( |4.5[ ) and the result 

(5.6) 


in Proposition 4.2 we get 


iiw = or). 

5.3. Estimating ( C w Here we will prove that 


^ h (£ u ~ Pu{K))^n,L 


T 2 (W) 


= o(h°°). 


(5.7) 


In view of the definition of the function in Proposition 4.2 we write, 


(5.8) 


+ e 


Xu{o)[£u,,x{ T 1t )\ + x(T 1 t)[£, u ,Xm{v)\ )^h,u , 


—1, 


where [•, ■] denotes the commutator. 
Then we have 


e * U 0 ,N,K/fi(£ u _ ^(h))^ = e ^W-^/ h 0 (h°°) 

+ e (^, 0 ^.n-^)/R^( a )[^ jX ( 7 -i T )]^ e ^/ft^ w + e (* u , w -*«)/fi 0LM(Mpp (!). (5.9) 
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Here the notation 0i / °o( suppx '.)(l) means that the function is supported on supp yt mid that it 

is uniformly bounded when h goes to 0. By Proposition 4.9, 0 < e^ w ' ai ' JV ' R— < 1 in w and 
for a € (0,1), 

&u},oj,N,h - < —(1 - Vl - 


in supp x'a,■ Now, (5.7) becomes a consequence of (5.9) and Proposition 4.2 thanks to support 
considerations. 


5.4. Proof of (5.3). Let us apply Proposition 4.7 with the following choices: u as in (5.5), 
z = l-iu(h) and 4> = ^u,u),N,h- We have 

ll e u llz, 2 (v w ) + w o ° e u Hl 2 (v w ) 

< Ch~‘\\e^^’ N ’ n/h (C UJ - Hcj(fr))u\\ L2( $^ + C'h _1 ||u|| Z/2( ^ ) . 

In light of (5.6) and (5.7), we deduce that, 




By Proposition 4.9, we have 

in 77 , and *.«)/» = 0(/r *) in L°° (A') 

In that way, we get the following estimate, 

l|e*"Nl mk) + l|S^®" /, «ll I .(*) = 0(S~), 

where 

K = IC x (0, T). 

We may rewrite this estimate in the form, 

\\e®“/ n u\\ L 2( K]L 2(Q^ + \\d a e^^ h u\\ L 2^ K . L 2^ T ^ = 0(h°°) , 


which in turn yields (5.3) in C(A; L 2 (0, T)) (cf. |6j Thm. 2; p. 302]). 
5.5. Proof of (5.4). 

Let (cf. (5.5)) ) 

v .— d T u — 

and 

w := d T v = d^u. 


We apply Proposition 4.7 to obtain, 


< Ch~ z \\e^ N ^ n {Cu - ^{h))v \| + C||i7| 


and 


w 


l\ e ^,*,N,K/h w \\ < cn- 3 ||e*«,ffl,.v,ft/ft(£ w _ MaJ (^)) w || + C]|i 

In light of the two identities 

Hw (h)) v — d T (^C*u — d T (^C-u Hu(fc)')'iph,Lu 

and 

(73^ A'oj( 7r))— d T (^C,u > 

we get by Proposition |4.2[ 

e ^ u / h {C ul - Huj(h)) v = 0 {h°°) and ~ Hu{K))w = 0(h °°), 


(5.10) 


(5.11) 


(5.12) 


(5.13) 

(5.14) 


(5.15) 
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By (5.6), we get further 


C u u = 0(h°°) in L 2 (%) 
Multiplying by u and integrating by parts yields, 


(\d T u\ 2 + h A a 2 \d a u\ 2 )adrda — f \u(a, 0)| 2 da = 0(h°°). 

J LU 


Note that the lowest eigenvalue of the operator — d 2 in L 2 (M + ) with boundary condition 
w'( 0) = — 2w(0) is equal to —4. Then we get 


_ 1 


and we deduce 


and / \d T u\ 2 drda — 2 / \u(a, 0)| 2 da > —4 / \u\ 2 drda, 

Ju JVu 


\d T u\ 2 adTdcr < 0(h°°) + 4|M|| 2 


In a similar way, using ( C w — ^(K^v = 0(h °°) (cf. (5.15|), we get, 


\d T v\ 2 adrda < 0(h °°) + 4||u| 


L 2 (V^) ’ 


Thus, we have the following two important estimates: 

IMI =0(h°°) and IHI =o(h°°). 

This and the estimates in ( ]5.15[ ) allow us to repeat the argument in Step 2 to obtain 

|| e ^,® I jv,n/fi( jCw _ ^(h))v\\ = 0(h°°) and \\e^ N ^ n {C u - ^(h))w\\ = 0(h°°). 


Now, (5.13) yields 


\\ e ^,w/n dT u\\ = o(h°°) and \\ e ^,w/^d 2 u\\ = 0(h°°) . 
As done in Step 3, the properties of the function N.h yield 

\\e^ /h d T u\\ = 0(h°°) and \\e^ /h d 2 u\\ = 0(h°°) in L 2 (K). 


Since |<h(J is bounded independently of h, then we infer from (5.11) 


e*“ /h d a u = 0[h°°) in L 2 (K). 

We recall that we have 

e*“/ R (£ w - ^{K))u = 0{h °°), in L 2 (K). 
Then, we use the estimates (5.16), ( 5.17| ) and 5.12 to deduce 

e^ /h dlu = 0{h°°) in L 2 (K). 

Since |<3?(J = 0(1), we get further 

= 0(h°°) in L 2 (K). 


(5.16) 

(5.17) 


This estimate and (5.11) yield the estimate in (5.4). 


6. Double mini-wells and interaction matrix 
In this section, we come back to the study of the double mini-wells operator C%. 
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6.1. Right and left operators. We introduce the operators corresponding to the left and right 
wells. Recall that we identify the boundary T = T with the interval (— L, L] and the orientation 
is chosen counter-clockwise. By Assumption |1.3| we know that 

{cr E (-L, L) : k(<j) = K max } = {—s 0 , s 0 } 

where sq E (0,L). We introduce the two (mini-)wells 

S£ = so and = — so , (6-1) 

and this is consistent with the counter-clockwise orientation of the boundary. 

Let us consider 77 such that 

0 < rj < ^ min (sg, L - ^ min ( - s r , -L + s r ) . 

We introduce the two intervals (in T) 

cog = {a E (— L, L] : |cr — s r | > 77 } and w r = {cr E (— L, L] : |cr — sg\ > 77 } . (6.2) 

We will apply the results in Section [4] with the interval u being ug or to r . The assumption on 77 
ensures that the left and right intervals oog and uo r have the same length 2 L — 2rj and are related 
by the simple transformation cr 1 —> —a. 

Let us introduce the two sets 


Wf = x (0, T), u r = w r x (0, T ), 
and the unitary transform U defined by 

Uf(a,r) = f (—cr, r). 

This transform goes from Li 2 {Cog) to L 2 (ca r ). Due to the symmetry assumption (cf. Assump¬ 
tion 1.3 1, we notice that 

Cg = UCrU' 1 . 

Thus these operators have the same spectrum and we may denote by 77 (h) their common lowest 
eigenvalue, i.e. 

fl(K) = flg(K) = Hr{K) • (6.3) 

The eigenfunctions of Cg may be deduced from the ones of C r . In particular, 

4>h,l = Ucf)n,T • (6.4) 


In Remark 4.4 we choose Co as 

Cog = {a : |cr — s r | > 277 } or Co x = {a : |cr — > 277 } . 

As a consequence, we get the two cut-off functions 

XI = Xuj e and %r = Xw r 

that are equal to 1 in Cog and Co r respectively. 

Proposition 5.1 yields, for every compact set K C Cog, 

= o{h °°), 

in C(K] L 2 (0, T)). Here 

= (V>, 4>hMw an d $v(cr) = $ r (-cr) = / yj K max - n(a) da . 


(6.5) 

( 6 . 6 ) 
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s = 0 



FIGURE 2. Illustration of the integral over a segment. In this case, we have 
I[a,b\ = fa alld f[c,a] = fc + f-L‘ 


6.2. Estimates of Agmon. We introduce the global weight 

d> = min (d> r , dv) , 

where 



We stress one more time that the integration over the segment [ 01 , 02 ] means the line integral 
along the boundary T from the point o\ to the point 02 hi the counterclockwise direction, see 
Figure 2. 

In particular, we have 

l*V)l a = «>(*)■ 


Let us define 

S u = <h r (s£), Sd = $e(s T ) and S = min(S u , Sd) • (6.7) 

Note that, on the “upper part”, - + d>£ = S u and on the “lower part” - + d>£ = Sd- In particular, 
we have 

d> r + dy > S . (6.8) 


The following proposition may be established by using the same estimates as in the proof of 


Proposition 4.7 


Proposition 6.1. For all a E (0,1), for all Co > 0, there exist positive constants ho, A, c, C 
such that, for all h E (0, ho), z E [—1 — At max ft 2 , —1 — At m «. Y ft 2 + Co ft 3 ], 116 H, 


f — z') u\\l 2 + Cft 3 ||u| 


L 2 {B(Ah 2 )) 


and 


ft 4 || 9 ( 7 ( e ^"^ u )||| 2 < Cft " 3 || e vT=H*/A ^ „|| 2 a + cfv 


3 "d' 2 


L 2 (B(Ah ?)) 


where B(g) = B T (g) U Be(g), where B r := B Ur , resp. Bj> := B Wi (cf Notation 4-6). 

6.3. Interaction matrix. 
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6.3.1. Preliminary considerations. 


Definition 6.2. Let us introduce the two quasimodes 

fh ,r = ~X.r ( t ) h,r i fh,C = Xlfih/ j 
that clearly belong to the domain of C} L . 


We use the following convenient notation. 

Notation 6.3. For M > 0, the notation 0(e~ M / h ) (introduced by Helffer-Sjostrand in [T8] ) 
stands for a quantity r(h,r]) defined on a set of the form (0, ho) X (0, 770 ) and satisfying the 
following: There exists a function 7 : (0, 00) —7 M such that lim^o l(jl) = 0, ant i f or all s > 0 
and 77 > 0, r(h,rj) = 0(e'- £ ^~' y ^~ M ^ h ) . The parameter 77 will measure the distance between \jio a 
and Co a , for a € {r, £}. 

The following lemma is the consequence of Agmon’s estimates and considerations on the 
supports. 


Lemma 6.4. For a € {r, £}, we let 

fh,a — {J-'h Kfr))fh,a — [£'TiX.a\(fh,a ■ 

Then, we have 

(i) r hta = 6(e~ s / h ), 

(h) (rn,a,fh,a) = 0(e~ 2S/h ) and (r h , a , fh,p) = 0(e~ s/h ) for a / /3, 

(hi) ( fh,a, h,a) = 1 + 0(e~ 2S/h ) and (f h , a , f n ,p) = 6(e~ s/n ) for ajC/3, 

(iv) If T = span {f h>1 , f ht g}, then dim T = 2. 

The following lemma states that the first two eigenvalues of £/> are close to p,(K) (the common 
first eigenvalue of the two mini-well operators) modulo 0(e~^/ h ). The proof is standard (see [4] 
or the presentation in 0). 

Lemma 6.5. Let us define Q = range l/ fi (£^) where In = (—00, —1 — K max h? + 27/i 3 ). Then we 
have 

(i) dist(sp (C h ),n(h)) = (5(e' s / n ), 

(ii) {(Cfr - fj,(h))u,u) > 'yh 3 \\u\\ 2 ,Vw € G L , 

(iii) dinrt/ = 2, 

(iv) sp(A) n I h c \p,(h) - 6(e- s / h ),n(h) + 0(e~ s l h )\. 


6.3.2. Interaction matrix. We want a more accurate description of the splitting between the first 
two eigenvalues of For that purpose, we will consider the restriction of JZn to the space Q 
generated by the first two eigenfunctions and we will exhibit an orthonormal basis of this space 
that allows us to compute asymptotically the eigenvalues of the corresponding 2x2 matrix. 

Let us introduce II the orthogonal projection on Q and g^ a = II f% a - As a consequence of 
Lemma 6.5 and of the spectral theorem, we get the following lemma. 

Lemma 6.6. We have in H 1 , 

fh,a ~ 9h,a = 6(e~ S/n ). 

From this lemma and Lemma |6.4[ we deduce the following. 


Lemma 6.7. Let us define the 2x2 matrix T by T a ,/3 = (fh,ai fh,p) f or a ft an d 0 otherwise. 
Then, we have 
(i) T = 6(e- s / H ), 

(h) = Id + T + d(e~ 2S / h ), 

(hi) (9h,on9h,p) = (fh,onfh,p) +0(e~ 2S/h ), 

(iv) ((gh,a,9h,p)) a ,p = Id + T + (5(e _2S / s ). 
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Let us now examine the so-called interaction matrix. The family (gn, a ) generates Q but is a 
priori not orthonormal. Thus we use the Gram-Schmidt matrix G = {{gn,at9hp))a,p and we let 
g = where g is the row vector (gn,e,gn,r)- The family g is now an orthonormal basis of Q. 

Let M be the (interaction) matrix of (the restriction of) in the basis g. 


Proposition 6.8. We have 

M = //(h) Id + W + d(e ~ 2S/h ), 

where W is defined by w a p = (r^ ai fh,p) if ol fi an d 0 otherwise. Moreover W is symmetric. 
In particular, the splitting between the first two eigenvalues of is given by 

A 2 (h) - Ai(h) = 2K, r (h)| + 6{e~ 2S / h ). 


6.4. Computation of the interaction. Now the problem is to estimate the interaction term 
Wi. r (h) given by 

w£ >r (h) = {(C H - n(h))f hj e,fh , r ) = ([£h,xe\<l>h,e,x r<0,r) ■ 

Let us recall that 

C r = — h 4 a~ 1 (d 0 -a~ 1 )d (T — a~ 1 {d T a)d T . 

Since xe does not depend on t, we get 

we, T (h) = -h A ([a~ 1 (d <T ^ 1 )d (7 ,xe\<l>h,e,XT<l>fv) % 

where 

xAA = Xr(-o-) • 

After the computation of the commutator and an integration by parts (with respect to a) to 
eliminate y”, we get 


= h i « X Xr Xe ({da4>h,r)4>h,i - (fhA^nA) dcrdr . 

Jv p 

Since y r = 1 in the support of y(, we get, 

W^ r (h) = h 4 / a^x'i ((<9<r<0r)<?V - (fn,r{da(fh,t)) dcrdr . 

Jv p 

Then, we integrate by parts and use the fact that is an eigenfunction of C a to get 

we, T (h) = wlfih) + , 

where 

Wg T (h) = h 4 [ a -1 r) - ^hA d ^h,e)} (0, r)dr , 

Jo 

wffih) = -h 4 f a -1 {(fn/idafinA - 4>hA d <r4>n,t)} (~L, T)dr. 

Jo 

Using Propositions 4.2 and 5.1 and the fact that 4>n,e( c r , r ) = <!>hA~cr, t), we write, 

exp 


ivlfih) = ^2h 5 / 2 |£ 0 ,£(0)| 2< f>^(0) + o(h 7/2 )) 
In the same way, we hnd 

w d e 


-S u 

h 


4AK) = (- h 5 / 2 toA-LAoAL){*AL) + &A-L)) + o(h 7 / 2 )) 


exp 


^Sd 

h 


this way, we get 


The computation of wf T is easy by using the expressions of £o,£ an d &e in Proposition 

XI.: st, 

- 2h5/2 0 1/2 L”® exp (- i 00 A + 0{nV2> \ exp (0) 


4.2 


In 


w e,r = 
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Note that t>(—L) = t)(L) by periodicity. In other words, we are saying that s = EL defines the 
same point on the boundary, see Figure 1. 

To compute wf T we use the two symmetry properties £o,e( a ) = £o,r( — cr) and <S>i(cr) = <h r (— a) 
and the expresions of £o,r and 4> r in Proposition 4.2 We obtain, 


w £,v = 


-2ft 5 / 2 0 1/2 v^M exp (- J° da ) + O(H^) 


exp 


-s u 

h 


By adding the expressions of w u and w d , we get an expression consistent with the one in (1.6). 
Recalling that h = /i 1 / 4 , we finish the proof of (|1.6|) by using Propositions 


3.3 


and 


6.8 


7. A Weyl formula 


This section is devoted to the proof of the following theorem (see [11] [233 f° r similar results 
for the Schrodinger operator with magnetic fields). 

7.1. Main result. For A £ R, we denote by 

N (£h, A) = Tr^/!/,, l(_ 0O ,A](^fe)) • 

the number of eigenvalues / x n (h ) of Ch below the energy level A. 


Theorem 7.1. Under Assumption \l . 1\ we have 

i. the Weyl estimate of the semiclassically negative eigenvalues: 

VA € (0,1), N (£ h ,-Ah) = ^^ + 0(1); 

/wo 7^2 

ii. the Weyl estimate of the low lying eigenvalues: 


VE e R, N (C h ,-h + Ehz 


1 


TTh* 


J r VW+ 


ds(x) 


(7.1) 


(7.2) 


tion 7.4 below). 


The proof of Theorem 7.1 relies on a comparison with an effective Hamiltonian (see Proposi- 


Remark 7.2. The counting of eigenvalues for the Robin problem appears (at least in the case 
of the disk) in the thesis of A. Stern [3?j in 1925 but note that the author (who refers to the 
book by Pockels |33j (written at the end of the nineteen-th century) is only counting the total 
number of negative eigenvalues. In this case, this is directly related to the counting function for 
the Dirichlet-to-Neumann operator. We refer to [13] for a recent survey on these questions. 


7.2. More about the Robin lD-Laplacian. This subsection contains one key element in the 


proof of Theorem 7.1 obtained through an additional analysis of the weighted operator in (2.7) 


and its groundstate. Note that the analysis of this operator is equivalent to that of the operator 


m 


- , s T \ 

B defined in (2.9). Recall that the operator R X B 1 is defined in the interval (0,T) and that 


U 

its ground state u B ‘ is given by the relation 

vSjP = (1 — Bt^vSjP 

where u^ is the groundstate of the operator in ( |2.7| ). Since u^ is normalized in L 2 (0,T), 

[ 8bu^ vS B ^ dr = 0 . (7.3) 

•Ai 

rjn 

For further use, we would like to estimate \\8bu b 1 1|z^ 2 ((0,T’),rfr)) uniformly with respect to B and 

T. 
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Lemma 7.3. There exist C > 0 and To > 0 such that for all T > Tq and B £ (—1/(3T), 1/(3T)) ; 

<r (7.4) 

(7.5) 


dnM [W B X ) <C, 
II^B^4 r} Hi 2 ((0.T),dr) < C . 


f - 

Proof. We recall that q B 1 is defined in (2.10) (the associated bilinear form is denoted in the 
same way). From the eigenvalue equation we get, for all p £ H l { 0, T), 


qPPP,p) ~ A i PPiv) = o. 

Then, we take the derivative with respect to B and we get 

q { B } (d B u { p ,p) - Ai (d B uP,p) = d B \i(u { p,tp) - d B q { B \u { B \ p) ■ 


(7.6) 


(7.7) 


We take p = vSp and the l.h.s. in ( |7.7| ) vanishes (use p = d B uP in ( 7.6| 1 ). Then, the r.h.s. 
vanishes and we deduce the Feynman-Hellmann formula 


d B Ai [n x B = 


'o 1,4(1 -BtY 


) up uPdr - ^\uP{0)\ 2 ■ (7.8) 


A T-uniform continuous Sobolev embedding (for T > 1) and Proposition 2.4 give 


\u B ^(0)| < C \\u B 4^(0,T) < C. 


(7.9) 


Therefore (7.4) holds thanks to (7.8), Proposition 2.5 and (7.9). 

- a* . r . — - - 

Now we take p = d B u x B 1 in (7.7) and we find with the same considerations (and (7.3)): 


qP(d B uP) - Ai (hiP^j \\d B uP\\ 2 L2 < C(\\d B uP\\ L 2 + \uP {0)\\d B uP {0)\) 

<c\\d B uP\\ L 2+c\d B uPm- 


With the spectral gap (see Lemmas 2.1 and |2.2| ) together with (7.3), we get, 

\\9buP\\ 2 L 2 <C + C\d B uP(0)\, 

and thus 

qP(d B uP) <C + C \d B uP(0)\ . 

From this we deduce 

II d B u B ^IIjjHo.t) < C + C || d B u B 4l 2 (o,t) + C I d B u B 4o)l 2 
and, by Sobolev embeddings, 

I d B u B 4o)| < C + C \\d B u B 4l 2 (o,t) ■ 

The estimate (7.5) follows from ( 7.10| ). 


(7.10) 


□ 


1 r j>i 

Note that ObUb 1 belongs to the form domain, but not to the domain of the operator. We have indeed 

(d B u { p (0) = d B (if 1 ) (°) 




2 
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7.3. Proof of the Weyl formulas. Thanks to the min-max principle and the usual Weyl 
formula in dimension one for the operator on the circle h? D 2 —K(a) (use a direct comparison with 


the case with constant potential for (i) and use for example [36] for the case (ii)), Theorem 7.1 


is a consequence of the following proposition which permits to localize the eigenvalues p, n of £h 
by comparison with effective Hamiltonians. 

Proposition 7.4. Under Assumption for eo £ (0,1), h > 0, we let 

A4 0 ,h = {n E N* : pL n (h) < -e 0 h} . 

There exist positive constants ho, C+, C— such that, for all h E (0, ho) and n E A f eo ,h> 

Tn{h) < Un(h) < n+(h ), (7.11) 

where n^(h) is the n-th eigenvalue of ’ defined by 


Cf' + = -h + (1 + C+h*)h 2 L% - K(a)h% + C+h 2 , 


and 


Cf~ = -/» + (!- C-hs)h?D% - k(ct)/0 - C-H 2 . 


7.4. Proof of Proposition 7.4 The proof will be done in three steps. 

7.4.1. Preliminary considerations. Thanks to the Agmon estimates established in Section [3] it is 


sufficient to work with C%. As suggested by the proof of Lemma 4.5, the spectral analysis of Cn 


may be done with the Born-Oppenheimer strategy. Let us recall the expression of the quadratic 


form Qfi, defined in (3.13), 
r L r T 


QhW = 

with T = DhT 1 . We let also 


f / a 2 ftf\d a ip\ 2 drda + 

-L Jo J-L 


\d r 'ip\ 2 d dr — \ip(a, 0)| 2 > da, 


u 


K,(cr),h 


= n {T} 


B 


with B = h* n(cr) = h 2 n(a). 

We introduce for a E [— L, L ) the Feshbach projection on the normalized groundstate of 
TJ K (a),ti> denoted by 

{'fji 'OK,[a),h)L 2 ((0,T)jadr)^ k((t), h • 

We also let 

= id - n CT 

and 

Rh{&) = \\daV K ^),h\\ 2 L^((0,T),ad T ) ' ( 7 - 12 ) 

The quantity is sometimes called “Born-Oppenheimer correction”. 

To be reduced to classical considerations, the main point is to control the effect of replacing a~ 2 
by 1. 

Lemma 7.5. We have, for all ^ E Dom(Q^), 

f d~ 2 \d a 'if\ 2 ddadr — f \d a ijj\ 2 ddadT 

lv T Jv t 

<c [ h 2 \f!^(a)\ 2 + hR h (a)\f^(a)\ 2 + h\d a U^i/j\ 2 dadT . 

Jv T 


with 


• (VK* 7 ? ')i V K,(a),h)L 2 ((0,T),adr) • 
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Proof. We write 


/ a 2 \d a if\ 2 da dr — \d a if\ 2 dadr 

lv T JV T 

< C f h^rld^iffdadr 

Jv T 

<2 C [ h 2 r (\d a n a if\ 2 + |5 CT n^r?/’| 2 ) dadr 


'Vt 


<c h I fJa) I + hR h (a)\f^{a)\ + dadr , 

Jv T 


where we used that 


r \v K (a),h\ 2 dr < C 


(7.13) 


(that is a consequence of Proposition 2.5) and that rh 2 may be estimated by Th 2 = Dh. □ 

Lemma 7.6. We have 

f T 

/ W(a-),?l < 5(Wre(<T) I 7). o dr 0(h ). 

J 0 

Proof. We notice from the normalization of v K ^fi that 


&0 / W(<x ) fiV n(a) jh® dr 0, 


so that 


2 f v. 
Jo 


r T 

K(a),hd*v K r a \ h ddT = / v K (a),H v K(a),h(d ( jd') dT , with d a d = ~Th 2 n'(a). 

Jo 


The conclusion follows from (7.13). 


□ 


7.4.2. Upper and lower bounds. Keeping these preliminaries in mind, the results below are con¬ 
sequences of almost the same computations as in |34.; Chapter 13] (see also |29j . [5] where a 


similar strategy is used). The first follows from a computation using Lemmas 7.5 and 7.6 


Lemma 7.7. There exist C > 0, ho > 0 siLch that, for all if E T>t and h E (0, ho), we have 
Qn( n,V0 < [ L + Ch 2 )\f^(a)\ 2 + {h\ 1 + Ch)Rn(a) + X,+ Ch 6 ) \U(a)\ 2 da . 


' —L 


The next lemma is slightly more delicate. 

Lemma 7.8. There exist C > 0, ho > 0 such that, for all if E T>t, e £ (0, ^) and h E (0,ho), 
we have 

QhW >J (l-e)(l-Ch 2 )h*\f^(a)\ 2 +{\i('H K{a)fi ) - C{e~ 1 h A R h {a) + e _1 h 8 + h 6 )} \f i) (a)\ 2 da 
+ r(l- e )(l- C 'n)n 4 ||5 CT n^||| 2(adT) +{A 2 (H K(CT) , ?l ) - Ce-WRrfa) - Ce~ l h 8 } \\U^if\\\, {ZdT) do. 

J — L 


Proof. First, we use Lemma 7.5 to get that 


QhW > Q n (if) - Chr / h 2 \f^(a)\ 2 + hR h (a)\f^,(a)\ 2 + h\d a Ilfif\ z dadT . 

Jv T 


(7.14) 


r L fT r L 

k4ia |2' 


Q^ PP (if) = / / K\d a if\adrda + 


l-L JO 


' —L W0 


\d T if\ 2 adT — \if(a, 0)| 2 > da. 


with 
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Then, we have the orthogonal decomposition 

f \d a ip\ 2 ad,T= I |n cr <9 cr '0| 2 a dr + f \U^: d a -i/j\ 2 a dr . 

Jo Jo Jo 

We have also the commutator identity 


(7.15) 


[da^a]^ L‘ 2 ((0,T),ad,T)'^K(cr),h ”1” (^5 ^/^(cr),/i)L 2 ((0,T),adT)^cr^/c(cr),/i 

- hi'(a)h 2 iff, r) rn K(tr)ih (r) dr^ , 

so that we get, by the Cauchy-Schwarz inequality, the estimate 

||[<9 ff ,n ff ]V>|| L 2 ((0JVdr ) < 2Rfi(a) 2 ||'i/’||L 2 (( 0 ,T),adr) + C7l 2 |[V’llL 2 ((0,T),adr) ■ (7-16) 

For all e E (0,1), we get, by using the classical inequality |a — b\ 2 > (1 — e)a 2 — e^ 2 b 2 and (7.15), 

f \d a ^\ 2 adr > (1 - e) { [ \d a Ii a ^\ 2 ad,T+[ \d a Ii^\ 2 a cLt\ - 2e^ x j |[<9 CT , n CT ]^| 2 adr . 

Jo U o Jo J Jo 

With ( |7.16 ), we get 


lo 


h 4 \d a 'i/j\ 2 adT > (1 — e)h A 


lo 


\da-Uuipf a dr + / |<9 (j n^V’l adr 

Jo 

,-1^4; 


- C'e- i (n 4 i? ft (a) + ) 11^11^2 ((0 , t ) ,adr) ■ (7-17) 

By computing and using Lemma 7.6 to deal with the double product, we have 

rT 

Z"ffc2\ | W ( \ 12 i / ~r> / /ot 2\ 


/0 


l&n^adr > (1 - Ch 2 )|/;(cr)| 2 + (i2 fi (<r) - Ch 2 )|/^(a)| 2 . 


(7.18) 


Moreover we have, by an orthogonal decomposition and the min-max principle, 

rT 

|2<j, 7_ !„/./■_ n\|2 \ \ /''!/ MX C_m2 , \ /^ / MItt-L„/.I|2 


/0 


1^1 adr- 0)1 > Ai(W /s( «r) ) ft)|/^(ff)| + X 2('H K (a),h)\\ Il a' l P\\L^(0,T),adr) ( 7 ' 19 ) 


The conclusion follows from (7.14), (7.17), (|7. 18), (|7. 19) and by integrating with respect to a. □ 


7.4.3. End of the proof of Proposition l.f We apply Lemmas 7.7 and 7.8 with e = h 2 . Then, we 
use Lemmas |2.1| and |2.2| and Proposition |2.4| to deduce that 

Ai(?4(c r),h) = -1 ~ h 2 n(a) + 0{h A ), 

and that there exist ho > 0 and C > 0 such that, for all h E (0, ho), 


eo 


A 2 {U K{a)fi )>-Ch>-^. 


Then we notice that Rfr{(?) (introduced in (7.12)) satisfies Rh(<j) = 0(h A ) thanks to Lemma 7.3 
and the relation B = K,(a)H 2 . The conclusion comes from the min-max principle (see (34] Chapter 


13]): the lower bounds in Theorem 7.1 follow from Lemma 7.7 and the upper bounds from Lemma 


Remark 7.9. One can see that Proposition 7-4 only requires that the boundary is C 2 and that its 
ciLrvature is Lipschitzian (that is an admissible boundary of order at least 3 in the sense of |32|). 
This residt matches with the one of |32| . Moreover, our effective Hamiltonians provide a uniform 
approximation valid for all the eigenvalues less than the energy level —e^h and not only for an 
h-independent number of low-lying eigenvalues. The underlying operator reduction follows from 
the general arguments often used in the Born-Oppenheimer framework. One can reasonably hope 
to extend the analysis to higher dimensional situations and improve the spectral approximations 
of |32| obtained in the case of admissible boundaries of order 3. 
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